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External field dependence of the correlation lengths in the 
three-dimensional 0(4) model * 
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We investigate numerically the transverse and longitudinal correlation lengths of the three-dimensional 0(4) 
model as a function of the external field H. In the low-temperature phase we verify explicitly the H~^^^- 
dependence of the transverse correlation length, which is expected due to the Goldstone modes of the model. On 
the critical line we find the universal amplitude ratio (,t/£,l = 1-99(1). From our data we derive the universal 
scaling function for the transverse correlation length. The //-dependencies of the correlation lengths in the high 
temperature phase are discussed and shown to be in accord with the scaling functions. 



1. INTRODUCTION 

In 0{N) spin models with > 1 two types 
of correlation lengths appear, corresponding to 
the transverse and longitudinal spin components. 
Like for the magnetization the behaviour of these 
correlation lengths in the critical region is de- 
scribed by asymptotic scaling functions and crit- 
ical exponents, which characterise the underly- 
ing universality class. In addition, there are pre- 
dictions |1I2| for the correlation lengths, which 
are related to the presence of massless Gold- 
stone modes. The measurement of the correla- 
tion lengths as functions of the external field H 
enables us to verify these predictions and to de- 
termine the critical parameters and scaling func- 
tions. We consider the three-dimensional 0(4) 
model because it is believed to belong to the same 
universality class as QCD with two degenerate 
light-quark flavours at its chiral transition in the 
continuum limit. The variant of the model, which 
we study here, is defined by 



pn = -J 
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where x and y are nearest neighbours on a hyper- 
cubic lattice with L points per direction, and Sx 
is a four-component unit vector at site x. From 



the lattice averages S*" and of the components 
parallel (longitudinal) and perpendicular (trans- 
verse) to H we find 



M = ( } , and ( ) = , 



(2) 



where M is the magnetization. Correspondingly, 
there are two susceptibilities 

^^'^ = V{{S\\^)-M^) , (3) 
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The connected correlation functions of the longi- 
tudinal and transverse spins are defined by 

Gi(x) = {Slsl)~M^ , Gt(x) = i(S^So^).(5) 

The large distance behaviour of Gl.t is governed 
by the exponential correlation lengths Cl.t, ex- 
cept for = 0, T < T^. 

The spontaneous breaking of the rotational 
symmetry for T < Tc gives rise to spin waves or 
massless Goldstone modes and leads to the fol- 
lowing divergencies for // ^ at all fixed T < Tc 

H-^'\ XT - H~\ - H-'/\ 



XL 



(6) 



*We thank for support by DFG Grant No. FOR 339/2-1. 



The prediction for comes from the relation 
XTj in accord with the mass interpretation 
— and ^ ^y^; the relation ~ xl 
(here, m — rria) does not hold for T < Tc PP. 
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2. NUMERICAL RESULTS 

Our simulations were done on lattices with pe- 
riodic boundary conditions and linear extensions 
L = 48, 72, 96 and 120. We used the Wolff single 
cluster algorithm and made 20000 measurements 
for each fixed (H, J)-pair. Between two measure- 
ments 100-3000 cluster updates were performed. 
In order to determine the correlation lengths we 
first calculated spin averages over planes and their 
respective correlation functions G(t). From the 
correlators at r and t + 1 we then derived an ef- 
fective correlation length ^e// (■'')■ When ^eff{T) 
reached a plateau inside its error bars we used the 
corresponding r-range to find ^ from a fit to 

G(t) - A[exp(-r/0 + cxp(-(L - t)/0] ■ (7) 

2.1. Critical behaviour 

In the thermodynamic limit {V — > oo) and 
close to Tc critical observables show power law 
behaviour in the reduced temperature t ~ {T — 
T,)/T, : For iJ = 

M = B{~tf , t < (8) 
XL = C+t-^ , t>0 (9) 
C = ^^t-" , t>0 (10) 

and at T = Tc for iJ > 

M = B^H^^^ , (11) 

= e2,T^""^ , '^c = ^/I3S . (12) 

Here, the temperature T is the inverse of the cou- 
pling J =1/T and Jc = 0.93590. From jSj we use 
the values B = 0.9916(5) and /3 = 0.380. We 
have fitted the data for M to Eq. ((TT|l and find 
B'' ^ 0.721(2) and S = 4.824(9). From /3, S 
and the hyperscaling relations all other exponents 
are fixed. In Fig. 1 we show the _ff-dependence 
of and at T^, and compare it to a fit to Eq. 
(|12|) . with the critical amplitudes 

= 0.838(1) , Cl - 0.421(2) , (13) 

which leads to the universal ratio = S.x/S.t — 
1.99(1), a result expected actually below Tc from 
a relation between Gl and Gt jH2j . At _ff = 
0, T > Tc, where = £,l = ^ and xt = XL 
we have made similar fits to our data and find 
the amplitudes ^+ = 0.466(2), G+ = 0.231(2). 
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Figure 1. The correlation lengths £,t and at 
Tc versus H from L = 48, 72 and 96-lattices. The 
lines are fits to Eq. jnj with the result ifT^ . 



From the measured amplitudes we obtain more 
universal amplitude ratios jl] 

i?x = G+{B'')-^B^-^ = 1.084(18) , (14) 
Qc = B^(^+)'^/C+ = 0.431(9) . (15) 
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Figure 2. The normalized transverse scaling func- 
tion gj{z) — S.rh'^" / gf {0) for various J-values. 
The line is the asymptotic form, Eq. ((T5)l . 

2.2. The scaling functions 

In the critical region the dependence of the ob- 
servables on the temperature and the external 
field is described by scaling functions 

M = h^l'faiz) , XL = h'/'-'f^{z)/Ho , (16) 
= h-'^^gjiz) , a = h-'^^g^iz) , (17) 
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where z = th'^/^^ and t = tB^/f^, h = H/Hq 
with Hq = {B'')~^ , are the normaUzed reduced 
temperature and field. The asymptotic behaviour 
of g^''^{z) for z ^ oo 0, t > 0), is 

9^i'^{^) ..^^^-^B^'"^-" ■ (18) 

In Fig. |2] we show the normahzed (and univer- 
sal) scaling function of the transverse correlation 
length. Whereas it was always possible to deter- 
mine ^T, we could only determine S^l above Tc, 
because below Tc a whole spectrum of states con- 
tributes to the longitudinal correlation functions. 
As it turned out, the shape of gj{z) is similar to 
foiz) and - at least for z > (T > Tc) - the 
shape of g^{z) resembles that of /^(z) (for ex- 
ample, the peak positions are the same). This 
similarity comes from the relations ^'^ ~ xt and 
Ci ^ XL- For more details see Ref. 5 . 

2.3. The i?-dependence at fixed T < Tc 

In Fig. O we show our results for S^t at fixed 
T = 1/J below Tc- The lines are fits to the form 
expected due to the presence of Goldstone modes 

£,T = xoR-^^^ + xi . (19) 

We see a clear verification of this prediction. The 
increase of the slope Xq = 0.481(2), 0.577(1) and 
0.637(2) for the three J-values with decreasing 
temperature is in accord with the scaling func- 
tion. By comparing the asymptotic behaviour of 
gj to the i7-dependence of we find 

^0 i-t)"^ , with q = -v + 135/2 = 0.1789 .(20) 
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Figure 3. in the low temperature phase at 
J 1.0, 0.975 and 0.95 vs. B'^l'^. The resuhs 
for J = 1.0(0.975) are shifted upwards by 10(5). 
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Figure 4. (filled) and (empty symbols) for 
J = 0.92 (o), 0.91 (O) and 0.90 (□). The resuhs 
for J = 0.92 (0.91) are shifted upwards by 5 (2.5). 



Indeed, the ratios between the different xq con- 
firm this formula. 

2.4. The iJ-dependence at fixed T >Tc 

Above Tci we have for all _ff > the inequal- 
ity > Cl, whereas at — the correlation 
lengths coincide. In Fig. 01 we show and t,L 
at T > Tc versus H. With increasing temper- 
ature the curves become flatter. The behaviour 
can again be explained from the asymptotic scal- 
ing function. Since ^ is an even function of h we 
must have 

^ = t-[go+git'^'''h^ + ...], (21) 

with constant gi. Yet, close to T^, where \z\ is 
small, but not /i, the expansion of g^ in powers of 
z leads to a different /i-dependence 

^ = h-'^[g^„+g'ith-'/^' + ...] . (22) 
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